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ABSTRACT 



Aims. ( Puppis, one of the closest and brightest massive stars, was the first early-type object observed by the current generation of 
X-ray observatories. These data provided some surprising results, confirming partly the theoretical predictions while simultaneously 
unveiling some problematic mismatches with expectations. In this series of papers, we perform a thorough study of ( Puppis in X-rays, 
using a decade of XMM-Newton observations. 

Methods. ^Puppis was observed 18 times by XMM-Newton, totaling lMs in exposure. This provides the highest-quality high- 
resolution X-ray spectrum of a massive star to date, as well as a perfect dataset for studying X-ray variability in an "archetype" 
object. 

Results. This first paper reports on the data reduction of this unique dataset and provides a few preliminary results. On the one hand, 
the analysis of EPIC low-resolution spectra shows the star to have a remarkably stable X-ray emission from one observation to the 
next. On the other hand, the fitting by a wind model of individual line profiles recorded by RGS confirms the wavelength dependence 
of the line morphology. 

Key words. X-rays: stars - Stars: early-type - Stars: individuals: f Puppis 



1. Introduction 
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With its very early spectral type (Q 4Infp, Wa lborn 
and a distance o f only 335 pc (Ivan Leeuwenl \10(J I ; 
Mafz Anellani z et all [2008), the star Naos, better known 
as £ Puppis (or HD66811), is one of the closest and brightest 
massive stars. It is therefore one of the most studied objects 
amongst the O-star population. However, despite the intense 
work, many open questions remain on its nature. 

Indeed, £ Puppis displays several intriguing properties. First, 
its visible spectrum shows clear signs of helium overabundance 
and chemical enric hment by CNO-processed material (e.g. 
Paul drach et al . 2001) as well as fast rotation (more than 200km 
s" 1 for v ship). iPennvl 119961: iHowarth et aljfl997l). Second, it 
is a k nown runaway (e.g. from Hipparcos data, Moffat et al. 
1998). These properties have led to speculations on its evo- 
lutionary status. On the one hand, the chemical enrichment 
and fast rotation could result from mass and angular mo- 
mentum exchange through Roche lobe overflow in a binary. 
£ Puppis could therefore have been the secondary component 
of such a system, the supernova explosion of its companion 
having ejected it from its bir th place a few millions years ago 
(van Rensbergen et al., 1996). On the other hand, £ Puppis dis- 
plays a similar Hipparcos pa rallax as stars of the Vela R2 associ- 
atiorQ ( Schaerer et al., 1997), and dynamical interactions within 



* Based on observations collected with XMM-Newton, an ESA 
Science Mission with instruments and contributions directly funded by 
ESA Member States and the USA (NASA). 
** Research Associate FRS-FNRS 

1 This conclusion was based on the original release of the Hipparcos 
catalog, hence the use of the 'old' distance of 430 pc in the Schaerer 



this a ssociation could have led to the ejection of the (single) O- 
star (Ivan R ensberge net all Il996l) . In this scenario, the chemi- 
cal enrichment of £ Puppis would be explained by the intense 
rotational mixing occurring in the fa st-rotating main-sequence 
progenitor (M evnet & MaedeR 120001) . In addition, f Puppis dis- 
plays doubl e-peaked emiss i on lines, suggested to arise in a ro- 
tating wind (IConti & LeepL 1 1974 iPetrenz & PulsL Il996l) . and a 
compression of the wind i n the equatorial plane was detected by 
lHarries & Howarthl d 19961) . 

Due to its brightness, £ Puppis was one of th e first mas- 
sive stars observed wi t h high -resolution in X-rays dKahn et all 
120011: ICassinelli et all 12001). At first, its X-ray lines appeared 
to match expectations as they did show the broad, blueward- 
skewed profiles expected for the wind embedded shock model 
dOwocki & Cohenl 120011) . However, the devil was in the de- 
tails. W hen quantitatively fitting the line profiles, Krame ret all 
(2003) found a much lower wind attenuation than expected on 
the basis of the mass-loss rate determined from optical and 
UV observations (see also [Oskinova et al.l 120061) . They also 
found that the typical optical depths r«, used in the wind- 
shock models, seemed independent of wavelength, which can 
only be explained by invoking porosity (Feldm eier et all l2003h 
lOskinova et all 120061). To improve the fitting of the X-ray line 
profiles. iLeutenegger et alj d2007l) included the effect of res- 
onance scattering: better fits were indeed obtained, without 
the need of a large reduction in the mass-loss rate, but they 
also showed that some unexplained dis crepancies remain. Re- 
analyzing the Chandra data of i Puppis. ICohen et alJ (1201 Oh ar- 



et al. paper, but the parallax similarity remains when using the new 
reduction of Van Leeuwen (and thus the 'new' distance of 335 pc). 
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gue in favor of a reduced mass-loss rate, without the need of any 
porosity as their new derivation of the optical depths implies an 
increase with wavelength, as expected from the bound-free ab- 
sorpti on opacity of the (cool) wind. Except for Leutenegger et al. 
(l2007h . all above studies relied on a single 68 ks Chandra obser- 
vation or a 57 ks XMM-Newton exposure taken in 2000. Both 
facilities have their advantages: while XMM-Newton globally 
has a higher sensitivity, Chandra has a lower background, and a 
higher spectral resolution and sensivity at short wavelengths for 
its grating spectra. Today, however, much more data are avail- 
able (see below). 

Considering the uncertainties in the line profile results and 
the lack of new variability studies, we decided to re-investigate 
^Puppis using the best dataset available at the present time: 18 
XMM-Newton exposures, corresponding to an exposure of lMs 
totaling >700ks of useful time (i.e. an improvement by an order 
of magnitude compared to most previous studies). This dataset 
thus provides the most detailed X-ray view of an O star to date. 
The results that we obtained will be presented in a series of pa- 
pers. The first one will present the data, their reduction, and a 
few first results; the second one will focus on the X-ray varia- 
tions of ^Puppis, using EPIC and RGS data; the last one will 
present a global analysis of the high-resolution X-ray spectrum, 
using the merged high-resolution data. 

This first paper is organized as follows. The dataset and its 
reduction are presented in Sect. 2, the spectral fits are presented 
in Sect. 3, the individual line profile fitting in Sect. 4, and the 
results are summarized and discussed in Sect. 5. 

2. XMM-Newton observations 

In the past decade of XMM-Newton observations, the star 
^Puppis was observed 18 times, mostly for calibration pur- 
poses. These datasets are excellent for studying the variability 
of ^Puppis since (1) the scheduled exposure times were often 
long (up to ~60ks) and (2) the observing dates probe weekly, 
monthly, and yearly timescales. Unfortunately, many observa- 
tions were affected by soft proton background flares, resulting 
in total exposure times reduced by about 30%. Total net expo- 
sure times for EPIC-MOS, EPIC-pn, and RGS amount to 579 ks, 
477 ks, and 751 ks. A summary of the observations is given in 
Table Q] The successive columns provide the dataset ID (obsID 
and revolution number); the date at mid-exposure (in the format 
dd/mm/yy + UT time and JD-2450000., calculated using the 
start/end times of the observations listed in the on-line 'observa- 
tion lokator'); the mode as well as the scheduled, performed, and 
effective (i.e. after cleaning flares) exposure time for both EPIC- 
MOS and EPIC-pn instruments; the scheduled, performed, and 
effective exposure time for RGS. An empty column indicates a 
discarded or unavailable dataset (see below). Note that the target 
was placed off-axis in two observations (5.95' off-axis in Rev. 
0731 and 1.1' off-axis in Rev. 0903). 

2.1. EPIC data 

The EPIC data were reduced with SAS vlO.0.0 using calibration 
files available on January 1 2011 and following the recommen- 
dations of the XMM-Newton tearrQ. Different modes (small win- 
dow, large window, full frame, timing) as well as different filters 
(thick, medium) and position angles were used for these observa- 
tions, resulting in a somewhat heterogeneous dataset. To ensure 

2 SAS threads, see 
http://xmm.esac.esa.int/sas/current/documentation/threads/ 



the most homogeneous analysis, hence a meaningful comparison 
between datasets, two decisions were taken. First, a few observa- 
tions were discarded: those when the instruments were not "on" 
(aka CAL CLOS ED), those totally affected by flares, those with 
the source appearing totally or mostly in a CCD gap, those with 
very short exposure times (<10ks), those taken in timing mode, 
and those using a unique combination of mode+filter. This trim- 
ming process results in a final dataset composed of 9 observa- 
tions taken with large window + thick filter and 6 observations 
taken with small window + thick filter for EPIC-MOS; 10 ob- 
servations taken with small window + thick filter, 5 observations 
taken with large window + medium filter mode and 4 observa- 
tions taken with large window + thick filter for EPIC-pn (see 
Table [TJ. Second, the extraction regions were chosen to be as 
constant as possible. A single source region was used for EPIC- 
MOS, whatever the mode, but two background regions were de- 
fined, one for each mode since it was not possible to extract the 
background on the same CCD chip as the source for the small 
window mode. For EPIC-pn, a single source region was used for 
all modes, as well as a single background region for the large 
window mode but four different background regions were nec- 
essary for the small window mode. Table [2] gives the position 
and shape of each of these regions. 

2.1.1. Pile-up 

^Puppis is rather bright in X-rays, with EPIC count rates of 
~2ctss~ 1 (for MOS) and ~6.5ctss _1 (for pn with thick filter). 
These count rates are at the pile-up limit from the XMM Users' 
handbook for the large window modes (i.e. 1.8 and 6cts s for 
MOS and pn, respectively) but well below the limits for the 
small window mode (i.e. 5 and 50ctss~' for MOS and pn, re- 
spectively). Some pile-up may thus affect our EPIC data taken 
in the large window mode. To see how severe the pile-up is, we 
performed several checks. First, we inspected the event files: no 
event with PATTERN-26-29 was found. The pile-up is thus 
moderate. Second, we run the SAS task epatplot: some small 
but significant deviation from the "no pile-up" configuration is 
detected for the large window data, especially those taken with 
the medium filter. 

To get rid of the pile-up, we could extract the data in an an- 
nulus centered on the source. This would require the annulus to 
be perfectly centered on the source, especially since the PSF is 
far from being symmetric. Getting the exact position of ( Puppis 
in the datasets is however an impossible task. Indeed, the detec- 
tion algorithm is also disturbed by pile-up, so that the position 
found in this way for £ Puppis is not accurate. For example, the 
pipeline-processed data from Rev. 1620 yields a position some 
2.5" away from the Hipparcos position of £ Puppis (this sepa- 
ration is the maximum found in the pipeline-processed data), 
while a dedicated run of the detection algorithm using only the 
PATTERN '=0 events gives a position that is only 1 .2" away. One 
would immediately think of using nearby X-ray sources associ- 
ated with well-known stars - sources which are less bright in X- 
rays (thus unaffected by pile-up) but still bright enough to get an 
accurate position in each dataset. However, such nearby sources 
do not exist in the neighbourhood of £ Puppis. A perfect cen- 
troiding of annular regions is thus impossible and using always 
the Hipparcos position for annular regions may alter the source's 
properties, hence the results of the variability study that we wish 
to perform. 
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Yael Naze et al.: A detailed X-ray investigation of f Puppis 



On the other hand, there is an alternative way in cases of 
mild pile-up: using only PATTERN-0 events. After such a fil- 
tering, we compared the PATTERN-0 spectra extracted in a cir- 
cular region with the spectra extracted in annular regions (hence 
free of pile-up) using the usual PATTERN=0-12 for MOS and 
PATT ERN=0—<\ for pn. This check was done by fitting simple 
2-temperatures models on LW+medium data from Rev. 0156, 
LW+thick data from Rev. 0731 and SW+thick data from Rev. 
1814 for EPIC-pn, and on LW+thick data from Rev. 0156 and 
SW+thick data from Rev. 1814 for EPIC-MOS. The compari- 
son is excellent for EPIC-MOS data: fluxes and count rates dif- 
fer by <1% and best-fit spectral parameters are within the errors. 
The remaining difference can be attributed to the slightly larger 
noise in the spectra extracted in annuli and from calibration dif- 
ferences. The comparison is less perfect for EPIC-pn, especially 
for the data taken with the medium filter: the flux differences 
reaches 6% in this case, and best-fit spectral parameters are at 
2-cr from each other. The pile-up thus still has a small influence 
on the EPIC-pn data taken in the large window mode, and those 
data should thus be considered with caution. 

2.1.2. The final files 

A final check was made on the data from Rev. 1620, which yields 
the most discrepant (2.5") position for £ Puppis, if we trust the 
pipeline processing. We first derived the position of £ Puppis 
from the PATTERN-0 data using the SAS task edetectchain, 
and then extracted the spectra using a circular region centered 
on that position. We compared these spectra to those extracted 
on PATTERN-0 data using a circular region centered on the 
Hipparcos position of £ Puppis. Both sets of spectra appear iden- 
tical in Xspec: the small centroiding errors have thus no impact 
on the spectra as long as a circular region is used. 

We therefore extracted lightcurves and spectra of £ Puppis in 
a circular region centered on the Hipparcos position of the target. 
We used only the PATTERN-0 event files. While this is not 
necessary for the small window mode, it ensures a homogeneous 
data reduction. Xspec vl2.6.0 was used to fit spectra, and our 
own software to analyze EPIC lightcurves. 

Note that spectra were grouped using the new SAS task 
specgroup. It enables to reduce the oversampling, which may 
"cause problems during spectral fitting because the spectral bins 
are then not completely independent" (excerpt from SAS 10.0.0, 
online documentation). We choose an oversampling factor of 5, 
ensuring that no spectral bin is narrower than 1/5 of the full 
width half maximum resolution at the central photon energy of 
the bin. Note that, while providing more statistically correct data, 
this process dramatically reduces the number of spectral bins. 
The data were also grouped to ensure that a minimum signal-to- 
noise of 3 was reached in each spectral bin of the background- 
corrected spectra. 

2.2. RGS data 

The RGS datasets were also reduced in a standard way with SAS 
vlO.0.0. Many new, important RGS features were modified in 
that version (e.g. the spectral binning in wavelength rather than 
in dispersion angle units). This ensures a better calibration of 
our datasets. It also solved the calibration problems (wavelength 
shift and reduced flux) found when using earlier versions of the 
SAS for the two observations where £ Puppis was placed off- 
axis. Note that the data were extracted using the proposal posi- 
tion of the source, which is the same in all observations but the 



first two (Revs. 0091 and 0156, shift of 0.0002° in both RA and 
DEC) - this small position shift has no impact on the derived 
RGS spectra. 

When detected, flares were discarded using rgs 

filteB The 

tasks rgsspectrum and rgsrmfgen then provided unbinned source 
and background spectra, as well as response matrices for each 
order (1,2) and each instrument (1,2). A final, combined spec- 
trum was also calculated using all 1 8 RGS datasets and the task 
rgscombine. The background files and matrix responses were at- 
tached to the source spectra using the new SAS task specgroup, 
which we also use to ensure an oversampling factor of maximum 
5 (see above). 

Fluxed spectra combining both RGS instruments and both 
orders were obtained using the task rgsfluxer. Note that a cor- 
rection for off-axis angles is applied to ensure that the fluxes 
are real photon fluxes and not simply recorded count rates (i.e. 
the arf response matrix is fully taken into account). The spec- 
tra of one revolution were sampled to get 1500 spectral bins, 
while the spectrum combining the 18 datasets was calculated to 
get 3000 spectral bins. This ensures an oversampling factor of 
about 3 and 6 for the former and latter cases, respectively. Two 
caveats should be noted. First, the rmf matrix is not fully taken 
into account by rgsfluxer and there is no instrumental width 
correction. This needs to be accounted for when modelling the 
spectra (see Paper III). Second, there are known small wave- 
length shifts in RGS spectra, apparently depending on the Sun 
aspect angle. However, no sign of such an effect is detected in 
our dataset when we use SAS vlO (they would appear as small 
spectral variations, and there are none, see Paper II). 

3. EPIC spectra 

With data of such high quality, the error bars on the spectra are 
very small, and it is therefore very difficult to get a formally 
acceptable fit. We thus avoided to try to get a perfect fit (i.e. 
X 1 ~ 1), which is actually impossible to get without going into 
unrealistic, overcomplicated models (e.g. 10 components fits 
with independent, free abundances) when all instruments agree. 
Rather, we have tried to get a fit as simple as possible which is 
at the same time as realistic as possible and as close as possible 

to t he spectral data. 

IZhekov & Palial J2007h showed that the high-resolution 
Chandra spectrum can be fitted by a shock model where the 
dominant temperature is 0. 1-0.2 keV Puppis was actually the 
star with the coolest dominant plasma), though some contribu- 
tion from plasma with 0.3-0.7 keV was needed to achieve a good 
fit. This modelling clearly shows that the plasma in £ Puppis 
is rather cool. We decided to use a simpler formalism, which 
does not try to reproduce shocks, but simply considers the ad- 
dition of optically-thin thermal plasma (without any assumption 
on their origins). In Xspec, the models using a distribution of 
thermal plasma (e.g. c6pvmkl) fail to provide a fit close to the 
data: we therefore had to fit the data using a sum of individ- 
ual optically-thin plasmas. As could be expected, one, two or 
three temperature fits do not provide fits close to the data, and 
similar conclusions are reached for four-temperature fits with a 
single absorption (which is unsurprising in view of the c6pvmkl 
result). Our model of choice is thus a four temperature model 
with individual absorptions: tbabs x 2 vphabs x vapec, where 
the first component represents the interstellar absorption, fixed 

3 Rev. 1071 has an increasingly high background towards the ob- 
servation's end, but no "discrete" flare. The whole observing time was 
therefore used. 
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Table 2. Regions used for extracting source and background data from EPIC instruments. 
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Notes. Here, 1 px=0.05" . The source position is from the Hipparcos catalog (cf. Simbad). The background regions Bkgd 2 to 4 were used for Revs. 
0636+0795+0980+1343, Revs. 0903+1071 + 1096, and Revs. 1620+1814+1983, respectively. In all other cases, the background region Bkgdl was 
used. 



to 8.9xl0 19 cm" 2 dDiplas & SavageUl994f) . and the abundances 
of the absorption and emission components are assumed to be 
equal. Adding a fifth thermal components (e.g. at 1 keV) does 
not significantly improve the fit, and we thus sticked to the deci- 
sion of using 4 components. 

After the best-fit temperatures and absorptions were found, 
we investigated the impact of using non-solar abundances. 
We began by keeping the helium abundance to 2 times solar 
(Rep olust et al.L 120041) . and let the nitrogen and oxygen abun- 
dances vary. We perfom such a simultaneous fit to all available 
EPIC spectra of each revolution. A few conclusions could be 
drawn from this trial. First, the pn spectra obtained in Large 
Window mode with the Medium filter (Revs 0156 and 0552, and 
second pn observation of Revs 535, 538 and 542) are clearly de- 
viant from other data, showing the impact of pile-up. Since all 
other data (pn or MOS, LW or SW modes with Thick filter) over- 
all agree with one another, we only discarded the pn data taken 
with Medium filter from further analysis. Second, the absorp- 
tions and temperatures of the fits do not vary much: we therefore 
fixed them to kT t of 0.09, 0.27, 0.56 and 2.18keV and N HJ of 
0.1, 0.1, 0.71, 0. xl0 22 cnr 2 . Finally, the origins of the high x 2 
can be better pinpointed. On the one hand, MOS and pn data do 
not always agree, especially at 0.4 keV (see Fig.Q]). This explains 
why some fits deviate from the mean behaviour: Revs 0731 and 
1071 only provided pn data, while Revs 0156 and 1164 con- 
sist only of MOS data. This difference between MOS and pn is 
probably due to remaining cross-calibration problems, not pile- 
up since f Puppis is far from the pile-up limit for data taken in 
SW mode with the Thick filter. On the other hand, all three in- 
struments (MOS1, MOS2, and pn) sometimes display 6x 2 of the 
same sign. This is often the case near strong lines, which points 
toward two explanations: (1) atomic parameters are imperfect 
(even for APEC - Fig.[TJshows two examples: at 1.24keV, there 
is flux in the data but not the model - a line is probably missing; 
at 1 .5keV, there is a line in the model but not the data) and (2) 
the asymetric line shape of £ Puppis influence even EPIC data. 
We cannot do much for the former, and the latter will actually 
be studied in detail in the third paper of this series - it is thus 
beyond the scope of this contribution. 

The next step is to free more abundances, but this could yield 
erratic and/or problematic and/or not better results. For exam- 
ple, the carbon abundance clearly gets unrealistic (whether freed 
last or first). Indeed, low-resolution, broad-band spectra such as 
those taken by EPIC yield few constraints on the carbon abun- 
dance and, as often happens with such data, the fitting proce- 
dure favors high values of carbon enrichment, whatever its ac- 
tual value. From previous studies, it is well known that nitro- 



Rev 1983: Nh and kT fixed 
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Fig. 1. The EPIC pn (top, in green) and MOS (bottom, in red and 
black) data of Rev. 1983 superimposed on the best-fit model. 



gen is overabundant in ^Puppis, and carbon may be subsolar 
(e.g. 0.35 times solar in Pauldrach et al. 2001 and 0.6 times so- 
lar in lOskinova et al.ll2007l) . However, reported abundances vary 
quite a lot in literature (se e other examples below) and, follow- 
ing |Zhekov&JPalla (2007), we decided to keep it to solar as for 
other non-constrained elements. For other elements, freeing the 
abundance may not improve the fitting quality or may not yield 
abundances significantly non-solar: this was the case of magne- 
sium, sulphur, iron and neon when released one after another. 
The silicon abundance stays close to solar but does improve the 
X 1 , it was thus allowed to vary together with He, N, and O. 

The best fit results are shown in Table [3] The parameter er- 
rors are taken from the raw fitting results of Xspec (i.e. these 
l-o" errors are "calculated from the second derivatives of the fit 
statistic with respect to the model parameters at the best-fit", and 
are indicative, see Xspec manual). No error is provided for the 
fluxes. We could indeed use the relative error on count rates, 
which would yield relative errors in the 0.1-0.3% range due to 
the sole Poisson noise. However, relative flux differences be- 
tween MOS and pn calibrations amount to about 1% in the soft 
energy band (where they are maximum), and the use of other 
similar models also yield 1 % relative errors. Fluxes should there- 
fore be considered as determined to 1/100, not 1/1000, uncer- 
tainties. 

Note that the listed parameters should not be over- 
interpretated: they simply represent a convenient way of well fit- 
ting the data, no more no less. The abundances, in particular, are 
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indicative, high-resolution data providing much more stringent 
constraints (see Paper III). It is for example interesting to note 
that freeing all abundances at the same time could lead to a better 
X 1 but to totally erratic and unrealistic abundances (carbon be- 
coming largely overabundant, nitrogen being solar). This being 
kept in mind, it is quite remarkable that our results, which should 
be considered as indicative only, agree rather well with previous 
abundance determinations. The helium abunda nce of ^P uppis 
was fo und to be 1.2 and 3 .4 tim es solaiQ by [Pauldrach et al.l 
(1200 11) and O skinova et al.l d2006l) . respectively, and we found 
an average value of ~2 (though this abundance is only weakly 
constrained in EPIC spectra). The nitrog en abundance was deter - 
mined to be 1 .7, 6, an d 8 ti mes solar bvlZhekov & Pallal (120071) . 
lOskinova etal] d2006l) . and lPauldrach et al.l(l2001l) . respectively, 
and we found again an average value of ~4. The oxygen abun- 
dance is the least constrained of all, with values of 0.75, 1.6 , 
and even 0.16 times solar reported by iPauldrach et ail (12001 1) . 
lOskinova etal] d2006l) . and lZhekov & Pallal (120071) . respectively. 
The latter agrees well with our value, but it is formally uncon- 
strained since its error was 0.23. O ur silicon determ i nation also 
agrees well with the value found bv lZhekov & Pallal (12007). 

A few general conclusions can be drawn. The dominant 
components are those with temperatures of 0.09, 0.27 keV, 
and 0.56 keV (providing 20%, 43%, and 35% of the to- 
tal flux, respect i vely), corroborating the conclusion found by 
Zhe kov & Pallal (1200%, with a different formalism, that cool 
plasma dominates in £ Puppis. The flux and spectral pa- 
rameters do not vary much, the largest variations are ob- 
tained when only one instrument is available (e.g. Rev. 0156, 
where the pn data are discarded because of the piled-up as- 
sociated with the Medium filter). Fluxes (and their associ- 
ated dispersions) in the total (0.3-4. keV), soft (0.3-0.6 keV), 
medium (0.6-1.2keV), hard (1.2-4.keV), and Berghofer's (0.9- 
2. keV) bandsO are 15.3±0.4, 4. 15+0. 14, 8.07+0.26, 3.06+0.09, 
5.67+0.17 xlO~ I2 ergcirT 2 s _I , respectively. Dispersions typi- 
cally amount to 3%, slightly larger than the typical 1 % error, and 
a shallow decreasing trend is detected, with a decrease of only 
~5% in flux since Rev. 0156. Such a decreasing trend is rem- 
iniscent of aging detector sensitivity problems, and we cannot 
exclude this possibility on the sole basis of the XMM-Newton 
dataset. 



4. Wind profiles 

ICohen et al.l (|2010) reported small variations of the line pro- 
files with wavelength, due to the energy-dependent opacity of 
the cool wind. However, the Chandra observation that they used 
was relatively short, hence subject to a much higher noise on 
the spectrum than for our data, and it is also rather insensitive to 
wavelengths >20A, where the effect is expected to be the largest. 
We therefore re-investigate here the issue using the combined 
RGS spectrum, which has a better signal-to-noise ratio and ex- 
tends beyond 20A. 

The left panel of Fig. [3] shows the observed Lyman a lines, 
in velocity space. For the figure, the lines were approximately 
continuum-subtracted and normalized to have a peak amplitude 
unity, to highlight their differences. Neighbouring lines can be 
seen for some of these Lyman a lines as bumps in the blue or 
red wings. Note that we do not show the N vn Lya line as it is 



4 As in Xspec: abundances are in number, relative to hydrogen, and 
relative to solar. 

5 For more details on the choice of energy bands, see Paper II. 
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Fig. 2. Spectral parameters as a function of Revolution number. 

blended with a N vi line. The variations with wavelength are ob- 
vious, as the peak velocity clearly appears less blueshifted for 
the short-wavelength lines. The comparison of width and skew- 
ness is more difficult by eye, as the RGS resolution broadens the 
short-wavelength lines in velocity space, blurring the trends. 

To quantify the wavelength variatio ns, we fit the lines with 
the same model^las lCohen et al.ld2010l) . to ensure homogeneity. 
Results are provided in Table [4] the first two columns identify 
the considered line, the next three columns define the line shape 
(characteristic continuum optical depth r», radius R* for the on- 
set of the X-ray emission, and the strength of the line), and the 
last two column provides details of the line ratios in the He-like 
fir triplets. 

Several things must be noted. First, in two cases (the He-like 
triplets of N and O), resonance scattering was needed to achieve 
a good fit. Second, the Lyman a line of nitrogen is blended with a 
line from N vi. We fit these two lines together, assuming that the 
line profile parameters are identical: keeping them independent 
yields unrealistic results (t* ~ 0) for the weak N vi line, and the 
lines are so blended that little independent information is avail- 
able, explaining the apparently strange results for the weakest 
line. The achieved fit of the nitrogen blend is far from perfect, 
however (x 1 ~ 2). Third, the iron line at 15A is not very well fit- 
ted (x 1 ~ 2), despite our efforts. It seems that line blends (there 
are numerous Fe xvin lines in the neighbourhood) affect the pro- 
file: though these lines are weak, the very low noise of our data 
reveals their impact, which was not obvious in the Chandra data. 
Finally, the fitting was done using a sin gle terminal velocity for 
the wind (2250km s -1 . lPuls et al.l l2006). a (3-1 exponent for the 
velocity law, and a power law of zero slope to represent the con- 
tinuum. 



6 Wind profile models for Xspec are available on 
http://heasarc.nasa.gov/xanadu/xspec/models/windprof.html 
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As it is always compatible with zero and as it yields no im- 
provement of the fits, no global line profile shift was applied, 
except for the iron lines near 17A where the improvement is sig- 
nificant with a shift of only -15+2 km s _1 . Wavelength shifts 
between instruments/orders were envisaged, but they yielded er- 
ratic results generally without improvement of the fit, sometimes 
large values (incompatible with wavelength calibration reports) 
for example for the iron line at 15 A, values compatible with zero 
within 3 sigma, and values for one instrument/order compatible 
with those of another instrument/order within 3 sigma. Cross- 
correlation also suggests null shifts between exposures and be- 
tween different instrument/order combinations of the same ex- 
posure. Due to the broadness of its lines, £ Puppis is indeed not 
the best source to find such shifts: to do so, the XMM-Newton 
calibration team uses sources with unresolved X-ray lines. 

As shown in the right panel of Fig. [3] the optical depth varies 
with w avelength. Our fitting thus confirms the ICohen et alj 
(2010) preliminary results by extending them to longer wave- 
lengths and by decreasing the noise on most lines (except the 
bluest ones). The wind of £ Puppis therefore is unlikely to be 
composed of clumps which are fully opaque at all wavelengths, 
as had been suggested for some porosity models. Indeed, op- 
tically thick clumps should produce a grey opacity that would 
solely be determined by the geometry of the clumps. In princi- 
ple, these optical depth variations may be used to constrain the 
m ass-loss rate, for a given star+wind model. This was attempted 
in lCohen et alJ d2010h . and they found that a reduced mass-loss 
rate with non-solar abundances provided the best-fit results. For 
comparison purposes, the same theoretical opacity is shown in 
Fig. [3] with a dotted line: below 20 A, it indeed provides a good 
fit, but for larger wavelengths, the agreement is less good be- 
cause of the nitrogen edge at 26A - an edge which appears very 
strong as the nitrogen abundance is enhanced for £ Puppis. This 
discrepancy may be solved by adapting the value of the abun- 
dances and mass-loss rate of £ Puppis, but this is beyond the 
scope of this paper. Clearly, additional modelling is needed to 
reproduce the full behaviour of £ Puppis in X-rays (and this will 
be done in Paper III). 

In contrast, the onset radius appears remarkably stable and 

confined in the 1-1.5/?, range, which is quite usual in massive 

■ 1 1 — - — — i 

stars (Giidel & Naze, 2009, and references therein). The sole ex- 
ceptions are the Nitrogen triplet and Lya lines, but it would not 
be surprising that, as these lines have emissivities that peak at 
rather low temperatures, it simply indicates a formation further 
out in the wind. 



5. Conclusion 

In the past decade, about lMs of data were obtained by XMM- 
Newton on £ Puppis. Of these, about 30% is strongly affected by 
flares, reducing the useful exposures to 579ks for EPIC-MOS, 
477ks for EPIC-pn, and 75 Iks for RGS. A variety of modes was 
used, the most reliable being the SW+Thick Filter mode; the use 
of the medium filter yielded piled-up data, while the data taken in 
LW+Thick Filter mode do not appear significantly different from 
those obtained with SW+Thick Filter mode. Attention was paid 
to this problem, notably by choosing similar extraction regions 
for all datasets and using only single events. 

Broad-band EPIC data taken with the Thick filter were an- 
alyzed using absorbed optically-thin thermal emissions. Four 
temperatures are needed to reproduce in a reasonable way the 
data. Note that the fits are not formally acceptable, but that (1) 
instruments do not always agree with one another and (2) the re- 



duced noise amplifies the limitations due to the imperfect atomic 
parameters and standard line profiles. Nevertheless, the EPIC 
spectra appear remarkably stable over the decade of observa- 
tions, with only 3% dispersions around the average fluxes. A de- 
tailed variability study, based on lightcurves, will be presented 
in Paper II. 

The combined high-resolution RGS spectrum confirms that 
the X-ray line profiles vary with wavelength. Fitting individual 
line profiles using a wind model yields similar onset radius for 
the X-ray emission, but wind continuum opacities depending 
on wavelength. This is simply due to the fact that the cool ab- 
sorbing clumps in the wind are not fully optically thick at all 
wavelengths, though further modelling is needed in order to ad- 
equately reproduce the opacity variations. 
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Velocity (km s -1 ) ^ 

Fig. 3. Left: Line profiles in velocity space, of the observed Lyman a lines. Right: Variations of the mean optical depth t„ with 
wavel ength. Resul t s from this work are shown with filled circles and hexagons - the latter for fits including resonance scattering, 
while ICohen et alJ d2010b results are displayed with empty triangles. The Cohen et al. theoretical absorption is shown by a dotted 
line. 

Table 4. Parameters of the wind profiles. 
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27 2.8 






7 6 93 


C vi-Lya 


33.7 


6 4 7 - 9 


1 - 1 i.o 
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Notes, tr refers to He-like triplets, and Ly to Lyman lines. Shifts and line profiles are identical for all lines of the doublets (Lyman a, Lyman 
p, Fexvn doublet at 17.1 A) and triplets; the line ratio in doublets is fixed to the ratio of maximum emissivities. For the Nvi triplet, resonance 
scattering parameters are t° > 56, < 0.28; for the O vn triplet, resonance scattering parameters are r° > 18, f} so i, = 1.4^. 
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